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ABSTRACT 

Since the radio-frequency emission from planets is expected to be strongly influenced by their interaction with the magnetic field and 
corona of the host star, the physics of this process can be effectively constrained by making sensitive measurements of the planetary 
radio emission. Up to now, however, numerous searches for radio emission from extrasolar planets at radio wavelengths have only 
yielded negative results. Here we report deep radio observations of the nearby Neptune-mass extrasolar transiting planet HAT-P-1 lb 
at 150 MHz, using the Giant Meterwave Radio Telescope (GMRT). On July 16, 2009, we detected a 3cr emission whose light curve 
is consistent with an eclipse when the planet passed behind the star. This emission is at a position 14" from the transiting exoplanet's 
coordinates; thus, with a synthetized beam of FWHM~16", the position uncertainty of this weak radio signal encompasses the location 
of HAT-P-1 1. We estimate a 5% false positive probability that the observed radio light curve mimics the planet's eclipse light curve. 
If the faint signature is indeed a radio eclipse event associated with the planet, then its flux would be 3.87 mjy±1.29 mjy at 150 MHz. 
However, our equally sensitive repeat observations of the system on November 17, 2010 did not detect a significant signal in the radio 
light curve near the same position. This lack of confirmation leaves us with the possibility of either a variable planetary emission, or a 
chance occurrence of a false positive signal in our first observation. Deeper observations are required to confirm this hint of 150 MHz 
radio emission from HAT-P-1 lb. 
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The importance of the physics behind star-planet interaction, 
in conjunction with theoretical predictions available for radio 
emission from a large number of extrasolar planets (GrieBmeier 
et al. 2007), has provided the impetus for making searches for 
decameter- and meter-wavelength radio emission from a few 
carefully selected extrasolar planets (Bastian et al. 2000; Ryabov 
et al. 2004; Winterhalter et al. 2005; George & Stevens 2007; 
Lazio & Farrell 2007; Smith et al. 2009; Lecavelier des Etangs 
et al. 2009, 2011; Lazio et al. 2010a, 2010b). Using the scal- 
ing laws (Zarka et al. 2001; Zarka 2007), radio emission from 
extrasolar planets orbiting close to their host stars could reach 
10 3 to 10 5 times the radio output from Jupiter and hence be de- 
tectable with the most sensitive existing radio telescopes. Even 
at large orbital distances (^1 AU), radio emissions from rapidly 
rotating planets orbiting stars with high X-UV brightness could 
be detectable (Nichols 201 1). Nonetheless, in the absence of de- 
tection, the theoretical estimates made for exoplanetary radio cy- 
clotron emission are based on a host of unknowns, such as the 
stellar wind, coronal density, and stellar and planetary magnetic 
fields. 

Up to now, various programs to search for radio emission 
have given negative results. The telescopes used include the 
UTR-2 (10 MHz, cr~1.6Jy), the VLA (74 MHz, cr~50mJy; 
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325 MHz, cr~0.6mJy; 1425 MHz, cr~15yuJy), and the Giant 
Meterwave Radio Telescope (GMRT; see review in Lazio et al. 
2009 and Zarka 201 1). Interferometric observations of high sen- 
sitivity and high resolution using the GMRT, which is a 30- 
km baseline Y-shaped array consisting of 30 steerable dishes 
of 45-metre diameter each providing an effective area around 
30 000 m 2 (Swamp 1990), seem very promising for this purpose. 

The planet HAT-P-1 lb is 4.3 Earth radii planet (Deming 
et al. 2011) orbiting a bright and metal-rich K4 dwarf with a 
4.89 day period on an eccentric orbit (e*0.2). This transiting 
planet located 38 pc away from us was discovered by Bakos et 
al. (2010) with ground-based photometry using the HATNet sur- 
vey (Bakos et al. 2004) and was confirmed by Dittmann et al. 
(2009). With a mass of about 26 Earth mass, this planet belongs 
to the so-called "hot-Neptune" category. The host star HAT-P-1 1 
is metal rich ([Fe/H]=+0.3 1+0.05). Its Kepler photometry has 
shown that the stellar surface has active spotted areas at fixed 
latitudes (Deming et al. 201 1; Sanchis-Ojeda & Winn 201 1). 

The planet's orbital axis is highly misaligned from the star's 
rotation axis, with a sky projected obliquity of A ~ 103° (Winn 
et al. 2010; Hirano et al. 2011). Nonetheless, the orbital plane 
is seen almost edge-on from the Earth and thus HAT-P-1 lb is 
eclipsed by the star for about two hours during each orbital cy- 
cle. The planet's eclipses can be used to discriminate any plane- 
tary emission from the stellar emission; the same technique has 
been used extensively for measuring planetary infrared emission 
using various space and ground-based telescopes (Deming et al. 
2005, Charbonneau et al. 2005, Sing et al. 2008). 
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2. Radio observations at 150 MHz and data analysis 

We observed the HAT-P-1 1 field with the GMRT before, during, 
and after the planetary eclipse, from July 16, 2009, 15h30m UT 
to July 17, lhOOm UT. The data consist of about 10 9 visibility 
measurements along the 435 baselines of the GMRT aperture- 
synthesis array. The integration time was 2 seconds and the cen- 
tral frequency was set at 150 MHz with a bandwidth of 16 MHz. 
The short integration time allows us to construct a radio light 
curve (see below), which is also sensitive to any large radio 
flares. 

Because the July 2009 run provided a ~3cr detection near 
the HAT-P-1 1 position and a radio light curve consistent with an 
eclipse of the planet HAT-P-llb (Sect. [3), we decided to carry 
out a second observation of the same field. This second run was 
scheduled on November 17, 2010, from 6h00m to 16h00m UT. 
Combining the two observation runs, the total data acquisition 
time on the target field is about 20 hours, including the total 
acquisition for the durations of the two eclipses (about 4 hours). 

The data reduction was done mainly using the AIPS++ pack- 
age. We used 3C48 as the primary flux density and bandpass 
calibrator, taking a flux density of 64.42 Jy at 148 MHz. The 
phase calibrator was the source 2038+513 which was observed 
repeatedly during the runs. 

While calibrating the data, bad data points were flagged at 
various stages. The data for antennas with relatively large er- 
rors in antenna-based gain solutions were examined and flagged 
over certain time ranges. Some baselines were flagged based 
on closure errors on the bandpass calibrator. Channel and time- 
based flagging of the data points corrupted by radio frequency 
interference (RFI) was done by applying a median filter with a 
6cr threshold. Residual errors above 5<x were also flagged after 
a few rounds of imaging and self-calibration. The final image 
was made after several rounds of phase self-calibration and one 
round of amplitude self-calibration where the data were normal- 
ized by the median gain found for the entire data. The final im- 
age was also corrected for the primary beam shape taken to be a 
Gaussian with FWHM of 173.83' at 150 MHz. 

After the final imaging stage, we generated light curves of 
several emission peaks detected around the HAT-P-11 position: 
one light curve of the peak source detected in July 2009 within a 
synthesized beam from HAT-P-1 1 (at such a low signal to noise 
ratio, the accuracy on the position of a detected source is ap- 
proximately the beam size), and the light curves of 14 emission 
peaks with similar S/N and detected around HAT-P-11 in both 
observations. These last light curves can be used for comparison 
purposes and validation of the statistics. To obtain these light 
curves, we calculated model visibilities using the sources de- 
tected in the entire field-of-view excluding a synthesized beam- 
wide region centered on the desired locations (ao, (5q); we then 
subtracted this model from the final calibrated visibility data. 
The residual visibility data (RVD) were then phase-centered on 
ao, #o and averaged for desired time bins to generate the light 
curves. 



3. Results 

3.1. Observations of July 16, 2009 

Figure [TJ shows the 150 MHz map of the field around HAT-P-1 1 
obtained from the observations on July 16, 2009. A faint source 
is detected at the coordinates a=19h50m51.0s, <5=+48°05'03" 
(J2000), which are 14" (about a beam width) north-east from the 
HAT-P-11 coordinates (a=19h50m50.247s, c5=+48 o 04'51.08"). 
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Fig.l. GMRT image of the HAT-P-11 field at 150MHz on 
July 16, 2009. The green ellipse in the lower-left corner shows 
the half-power beam width. The contour levels given at the 
top of the image are in units of Jy beam 1 and are defined 
as mean+rmsx(n), where values of n are given in the brack- 
ets. Thus, the first, second, and third contours correspond to 
1.5x1.05 mJy, 2.1x1.05 mJy, 3.0x1.05 mJy, etc. Negative con- 
tours appear as dashed lines. 



Figure |2] shows the 150 MHz light curve derived at the source 
position. 

Over the full integration time, including the two-hour 
eclipse, the source has a mean flux of 2.45+1.09mJy (2.2cr de- 
tection). If we consider only the time outside the eclipse (Fig. [2), 
a mean flux of 3.87+1.29mJy (3cr) is measured. During the 
eclipse of the planet behind the host star, no emission was de- 
tected and the mean flux was found to be -0.92+1. 94 mJy, con- 
sistent with zero emission for this time interval. The light curve 
is thus in agreement with the expected light curve for a planetary 
radio emission, eclipsed when the planet is known to be behind 
its host star. Thus, a straightforward, albeit heuristic, interpre- 
tation of the light curve would be a 3<x detection of the planet 
HAT-P-1 lb at 150MHz. 

The position of the source exhibiting the eclipse light curve 
has been measured to be 1 4" away from HAT-P- 1 1 . However, for 
such a faint source the positional accuracy is similar to the beam 
size of 16". More precisely, from our experience with GMRT 
maps at 150 MHz, we derived an empirical expression for the po- 
sitional accuracy on GMRT maps; at 90% confidence we found 
an accuracy on the position given by Beamsizex(l / S NR+ 1/10), 
which is in agreement with the same relationship found on the 
VLA First Survey Catalog which reports^ a 90% confidence po- 
sitional accuracy of Beamsize x (l/SNR + 1 /20). Applying this 
relationship to our measurement of the emission peak at 2.2<x 
and taking into account the 4" pixel size, we find that the above 
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Fig. 2. Time series of the 150 MHz flux density measured on July 
16, 2009, in the direction of the radio source near HAT-P- 1 1 . The 
measurements have been rebinned to 36 minutes. Triangles and 
squares correspond to the RR and LL polarizations, respectively. 
The two vertical long-dashed lines indicate the beginning and 
the end of the planet's eclipse behind the host star. The dashed 
horizontal lines show the box-shaped eclipse light curve fitted to 
the data points. 
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Fig. 3. Same as in Fig. [2] in the direction of HAT-P- 11 on 
November 17, 2010. The mean of the flux density measurements 
is found to be 0.00+1.19mJy. Thus, no emission is detected to- 
wards or near HAT-P- 1 1 in these observations. 



measurement of the radio position can be consistent with the 
HAT-P- 11 coordinates. 

To assess the reliability of this detection, we also estimated 
the false-positive detection probability, i.e. , the probability of ob- 
serving a light curve of this shape for a constant radio emission 
with no real eclipse, but only an eclipse mimicked by the ran- 
dom noise along the radio light curve. This is calculated by es- 
timating the probability that the mean flux during the eclipse is 
lower than zero, as observed. For this calculation we used the 
present estimate of the source flux, and the error bars for the 



Fig. 4. Time series of the 150 MHz flux density measured in 
the direction of one comparison source at a=19h50m32.2s, 
<5=+48°04'22" (J2000), at 183"from HAT-P-11. The dashed 
horizontal lines show the mean values of flux density at the two 
observation epochs. 
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Fig. 5. Same as Fig. [4] but for the observed peak at 
a=19h50m45.8s, (5=+48°04'10" (J2000), at 61"from HAT-P- 
11. 



times outside and during the eclipse. For the case of no real flux 
variation, our analysis shows that the probability of the mean 
flux during the eclipse being below zero is only 4.8%. We con- 
clude a 5% false positive probability of eclipse detection in the 
observed 150 MHz light curve at the position coincident (within 
the errors) with the planet HAT-P- 1 lb. 

We also computed the x 2 of two different models fitting the 
light curve data with a 339-second sampling time providing a 
total of 190 measurements. In the first single-parameter model, 
we assume a constant flux level for the entire light curve, corre- 
sponding to the hypothesis that the radio source is not the planet 
and therefore the flux density is constant in time. We find a x 2 
of 186.6 for 189 degrees of freedom. In the second model, we 
assume that the source is a planetary emission and the flux den- 
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Fig. 6. Same as Fig. S but for or=19h50m56.9s, <5=+48°05'46" 
(J2000), at 87"from HAT-P-1 1. 
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Fig. 7. Same as Fig.S but for cr=19h50m58.9s, (5=+48°02'26" 
(J2000), at 169"from HAT-P-1 1. 



sity during the eclipse is fitted by a second parameter, larger than 
or equal to zero. With this two-parameter model we find a x 2 of 
182.8 for 188 degrees of freedom. Using the F-test, we conclude 
that the decrease of x 2 by 3.8 favors the second model with a 
95.0% confidence level, which is consistent with the conclusion 
based on the false positive probability calculated above. 

To make the above analysis more rigorous, we also derived 
the radio light curve for 14 sources seen within a few arc minutes 
around HAT-P-1 1 and detected at both observation epochs. The 
measured fluxes of these comparison sources are in the range ~1 
to 3 mJy. Four examples of light curves are given in Figs. |4]to [7] 
We found that the false positive probabilities computed for these 
directions are uniformly distributed between and 1, validating 
our approach. Thus, to a 95% confidence limit, the detected sig- 
nal is consistent with a 150 MHz radio emission associated with 
the extrasolar planet HAT-P-1 lb. 



3.2. Observations of November 17, 2010 

In the full-synthesis GMRT observations of November 17, 2010, 
we did not detect any emission toward HAT-P-1 1 or nearby di- 
rections within a synthesized beam. On this occasion, the near- 
est detected radio source lies at 27.5" from HAT-P-1 1, and has a 
flux of 1.78+1.01 mJy. The light curve of this source at 150 MHz 
does not show any eclipse-like signature. Therefore, we are led 
to conclude that in these observations of HAT-P-11 we did not 
detect the radio source and its associated eclipse event we had 
seen on July 16, 2009. 

3.3. False positive or variability? 

There are two plausible explanations for the discrepancy be- 
tween the results of the two epochs. First, a false positive detec- 
tion at the first epoch cannot be completely ruled out given the 
estimated false positive probability of 5% (Sect. [XT). Moreover, 
the source is close to the detection limit. So a false detection due 
to noise cannot be excluded at present. 

The second interpretation invokes the general variability of 
planetary radio emission at low frequencies. The variable nature 
of the radio emission is indeed predicted by models including the 
observed variability of the stellar magnetic field. For instance, 
in the case of HD 189733b Fares et al. (2010) showed that the 
observed variations in the position of the so-called "source sur- 
face" (the stellar surface beyond which the magnetic field be- 
comes purely radial) produce variations in the expected radio 
flux. Thus, in the case of HAT-P-1 1, if the discrepancy between 
the results of July 2009 and November 2010 is due to variabil- 
ity of the radio emission, we may have actually detected spo- 
radic emission from the extrasolar planet at 1 50 MHz on July 1 6, 
2009. However, the variable nature of the emission precludes a 
claim of detection; we are left with a hint, and we must wait for 
further, possibly even more sensitive observations. 

4. Ephemeris uncertainties 

The characteristics and ephemeris of the planet's orbit are 
needed to compute the eclipse light curve. Here we used the or- 
bit's eccentricity e and the longitude of periastron co provided in 
the discovery paper by Bakos et al. (2010): e cos co-Q. 20 1+0.049 
and e sin <y=0.05 1+0.092. The more recent observations of 
transits of HAT-P-llb using the Kepler satellite provide new 
ephemeris for the transit central time To and orbital period 
P (Deming et al. 2011, Sanchis-Ojeda& Winn 2011). We 
used the following values given in Deming et al. (2011): 
7o=2454605.89155 (BJD) and P=4.8878018 days. The differ- 
ence between the eclipse time calculated using these new values 
and the time calculated using the To and P given in Bakos et al. 
is not more than 1.5 minutes. 

More importantly, the new measurements using the Kepler 
photometric data yield a smaller impact parameter b\ = 
acosi/Rp, where a is the semi-major axis, i the orbit's incli- 
nation, and R p the planet's radius. A smaller impact parameter 
leads to a longer eclipse duration. Using the new values, we 
found an eclipse duration of 2 hours and 27 minutes, which is 
38 minutes longer than the previous estimates using the Bakos 
et al. (2010) impact parameter. 

By far, the largest uncertainty in computing the eclipse time 
arises from the uncertainty in the orbital eccentricity and longi- 
tude of the periastron, the uncertainty reaching up to +3h30m 
in the extreme cases. Unfortunately, the eclipse amplitude in 
the infrared may be too faint to be detected using the Spitzer 
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Space Telescope (Deming et al. 201 1). Only new radial velocity 
measurements in the coming years may improve the accuracy of 
the orbital eccentricity. Here, we used the best published values 
given by Bakos et al. (2010) and used the eclipse phase 0=0.628, 
where <p-Q corresponds to the central transit of the planet in front 
of the host star. 

5. Discussion 

5.1. Radio emission frequency 

The principal mechanism put forward for planetary radio emis- 
sion is electron-cyclotron maser radiation. It occurs at the local 
gyrofrequency f g = 2.8 (b p /1 g) MHz, where B p is the planet's 
magnetic field. The detection at a frequency of 150 MHz would 
then correspond to a planetary magnetic field of around 50 G. 

We recall that the Jovian magnetic field estimated from the 
radio frequency cut-off in its cyclotron emission is around 14 G. 
Thus, the present tentative radio detection requires a magnetic 
field for HAT-P-1 1 b (~50 G) that is of the same order of mag- 
nitude, albeit slightly stronger than that of Jupiter. We have no 
clues to the possible origin of a magnetic field of this strength, 
although we can emphasize here the particularities of this planet 
which has a significantly higher density than Jupiter, and orbits 
a metal-rich star ([Fe/H]=+0.3 1+0.05). The planet HAT-P-llb 
is consequently expected to have a high abundance of heavy el- 
ements and a large core mass (Guillot et al. 2006). Taken all 
together, its peculiarities (eccentric and highly inclined orbit, 
metal-rich host star, short orbital period, etc.) make this planet 
different enough from any planet in our own solar system that 
it is difficult to make a reasonable guess of its magnetic field 
strength. Hence the requirement of a magnetic field 3 to 4 times 
Jupiter's seems quite plausible. 

5.2. Emission flux 

Several mechanisms for exoplanet radio emissions have been 
considered by GrieBmeier et al. (2007): stellar wind (SW)- 
magnetosphere magnetic interaction, SW-magnetosphere kinetic 
interaction, and coronal mass ejection-magnetosphere kinetic in- 
teraction. Application to HAT-P-llb of the scaling law extrap- 
olations of GrieBmeier et al. (2007) lead to predicted radio flux 
densities of 14 mJy for the SW-magnetosphere magnetic interac- 
tion, 1 mJy for the coronal mass ejection (CME)-magnetosphere 
kinetic interaction, and <<1 mJy for the SW-magnetosphere ki- 
netic interaction (note that these are mere orders of magnitude). 
The first two values are of the same order of magnitude as the 
~4 mJy flux density tentatively detected here. 

Alternatively, if the stellar magnetic field is strong enough, 
the planet may induce electron precipitations and thus radio 
emission in the star's magnetic field (Zarka et al. 2001; Zarka 
2007). This may be the case for HAT-P-1 1 b whose semi-major 
axis is only 0.053 au. The stellar magnetosphere could even ex- 
tend beyond the planets orbit. In the framework of the Jardine 
& Cameron (2008) model in which the radio emission is caused 
by interaction of the planetary magnetosphere with the stellar 
corona in which it is embedded, we can derive typical values for 
the resulting radio signal from HAT-P-1 1 b. Because HAT-P-1 1 
is a slowly rotating star (Bakos et al. 2010) with a rotation pe- 
riod of 30.5+1, 1 days (Sanchis-Ojed & Winn 2011), the planet's 
velocity relative to the star's magnetosphere is v p - v*=99 km s _1 
(i.e., three times higher than that of HD 189733b). Taking a 10% 
efficiency for conversion of the power of the accelerated elec- 
trons into radio emission within an emission beam solid angle 



of 1 .6 steradian, a stellar magnetic field Z?,=40 G as measured 
for the K-type star HD 189733 (Moutou et al. 2007), and a stel- 
lar density at the base of the corona equal to the solar value, 
we obtain a predicted radio flux in this model of the order of 
13 mJy from HAT-P-1 1 b, again similar to the tentative detection 
reported here. 

We conclude that a radio emission of the hot-Neptune HAT- 
P-1 1 b at a level of a few milli-Jansky level is consistent with ex- 
isting theoretical predictions, in particular the predictions based 
on the cyclotron maser instability (CMI) of energetic electrons 
accelerated through the interaction of the planetary magnetic 
field with the stellar magnetic field, or coronal mass ejections. 

6. Conclusion 

Although we are not able to draw a definitive conclusion on 
150 MHz radio emission from HAT-P-1 1 b, at the very least the 
hint of radio detection presented here identifies HAT-P- 1 1 b as 
a prime candidate for many follow-up observations in the near 
future. The priority is to try to confirm the present tentative de- 
tection, via re-observation with GMRT at 150 MHz and/or new 
observations with LOFAR in the 30-250 MHz range (providing a 
1-100 mJy sensitivity depending on the observation parameters 
used; van Haarlem et al., submitted, 2013) and UTR-2 in the 
10-30 MHz range (providing a sensitivity of ~100 mJy; Ryabov 
et al. 2004). The observed spectral range can also be extended 
toward shorter wavelengths using GMRT, with even higher sen- 
sitivities of ~1 mJy at 240 MHz and -50 /iJy at 614 MHz (e.g., 
Lecavelier des Etangs et al. 2009). Observations over a broad 
range of frequencies will, with any luck, allow confirmation of 
the existence of the emission and better constraints on the plane- 
tary magnetic field strength and determination of the radio spec- 
tral index. 

The future observations should also be distributed at multi- 
ple epochs and at different orbital phases of the planet to charac- 
terize the suspected variability of the radio emission in terms 
of duty cycle, as well as the radio emission as a function of 
the star-planet angle and the emission directivity (CMI produces 
narrowly beamed radio emission). In this context, Hess & Zarka 
(2011) have analyzed all the observables that could be derived 
from a broadband dynamic spectrum with a sufficient signal-to- 
noise ratio. 

In the long term, it should be possible to investigate if the 
radio emission from extrasolar planetary systems correlates with 
the stellar spot activity. In this context, HAT-P-1 1 represents the 
target of choice because the Kepler mission has recently demon- 
strated the possibility of following the spot activity and perhaps 
even of drawing the butterfly diagram of surface spots for this 
star (Deming et al. 201 1, Sanchis-Ojeda& Winn 201 1). 

Radio observations of different types of exoplanets are also 
crucial for comparative exoplanetology. For instance, observa- 
tions of the other transiting hot-Neptune GJ436 b (Gillon et al. 
2007) will allow comparison with a planet similar in mass but 
orbiting an M-dwarf star which is half the mass and radius of 
HAT-P-1 1 and known to have a strong Far-UV and Lyman-a out- 
put (Ehrenreich et al. 2010). With a 6.55 Earth mass, the nearby 
exoplanet GJ1214b should also provide an opportunity to un- 
dertake search for radio emission from a water-rich lower mass 
super-Earth planet (Charbonneau et al. 2009; Desert et al. 201 1). 
In the near future, more sensitive observations at even lower fre- 
quencies will become feasible with the LOFAR observatory, pro- 
viding a great boost to exoplanetary research. 
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